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1. Introduction 

Recent years have witnessed a lot of interest in the non-perturbative understanding of the 
symmetry breaking sector of the standard model. An important question is whether the 
scalar field self-coupling and the Yukawa coupling are trivial as suggested by the signs of 
the perturbative (3 functions. If the model is influenced by the Gaussian fixed point for 
all allowed values of the bare couplings, one can derive non-perturbative upper bounds on 
the renormalized Yukawa coupling yn and quartic self-coupling A/j from the constraint that 
the cut-off must be sufficiently larger than the masses of the particles. Using the relations 
niH = V^^rVr and mp = HrVr this is equivalent to finding upper bounds on the Higgs mass 
mH and the heavy fermion mass irtp {vr ~ 246 GeV is the scalar field expectation value). 

It is reasonable to address these questions first in simplified lattice fermion-Higgs models 
without gauge fields. For recent reviews on such models we refer the reader to ref. |jl|. 

In a model with naive fermions the number of mass degenerate doublets is as large as 
16. With the usual staggered fermions one can reduce this number of doublets to four and 
using the mirror fermion model with Wilson fermions of ref. it appears to be possible 
to reduce the number of doublets to one, while the 30 doublers acquire masses of order of 
the cut-off 0. For numerical investigations of these models with the hybrid Monte Carlo 
algorithm (HMCA) it is necessary to double these numbers of fermion doublets. 

In this letter we shall investigate a new proposal 0] for a fermion-Higgs model which is 
based on the reduced staggered fermion formalism and which describes one SU(2) doublet 
in the scaling region (two doublets when investigating the model with the HMCA). In this 
paper we shall introduce the model and present some preliminary results. A detailed account 
of our investigation will appear elsewhere |^ . 

2. The model and its symmetries 

The usual euclidean staggered fermions on a four-dimensional hypercubic lattice describe 
four flavors in the scaling region. By using the 'reduced' staggered formalism [|, |^ this 
number can be halved to two. If two of these staggered fields are placed in an SU(2) doublet 
one would get a model with two degenerate isospin doublets in the scaling region. In this 
letter, however, we follow another strategy. In order to use the fermion field as efficiently as 
possible we want to couple the two staggered fiavors to the Higgs field, leading to a model 
with only one doublet in the scaling region. 

To show how the Higgs field can be coupled to the staggered fermion fiavors while pre- 
serving as much symmetry as possible, we first write down the target action that we want 
to reproduce with the lattice model, 

Sf = - [ d'xii^rd.i^ + yi^L^^R + V^«0Vl)] (1) 



The Dirac doublet ip interacts with the scalar field (p G SU{2) which contains the 0(4) 
components (y9^ of the Higgs field in the familiar way, cj) = i J2^=i ^j^j + ^aX where Tj are the 
usual Pauli matrices. 

In order to reproduce this action with staggered fermions, we need to exhibit their spin- 



flavor structure and it is convenient to introduce the 4x4 matrix fields P] 



b b 

Here x ^^^^ X ^-re the usual (one-component) staggered fields, the sums in eq. (0) are over 
the 16 corners of the unit hypercube, 6^ = 0, 1 and 7^+'' is a short hand notation for the 
product 71^'*' ^72^"*" ^73^''' ^74"''' "• Since \l/ contains 16 times as many degrees of freedom as 
X, not all components "^""^ are independent. However, the components of the low momentum 
modes "^{p) with — 7r/2 < p^ < 7r/2 are independent [Q. Spin-flavor transformations on x 
correspond to discrete transformations on \1' from the left or right. A shift transformation 
Xx -^ CpxXx+p, to be interpreted as a discrete flavor transformation 0, |^, translates into 
"^x -^ "^x+plp and a 'spin' transformation Xx -^ VpxXx+p translates into "^x -^ 7p^x+p- The 
sign factors rj^x and C^a; are defined as r]^x = (— l)^i+"+^>-i and (^x = (— l)^>+i+"+^4_ 
The kinetic part of the free staggered fermion action 

Sk = -lJ2Vf^x iXxXx+p - Xx+pXx) (3) 

XjJ, 

can now be written as 

Sk = -Y1 iTr(^,7^^,+^ - ^,+^7^^,), (4) 

which reduces in the classical continuum limit to the gradient term in eq. (|l|), however with 
two doublets. 

In the reduced staggered formalism the field x is restricted to the odd sites Xx —>■ ^(^ — 
^x)Xx and Xx to the even sites Xx ~^ ^(^ + ^a;)Xa;5 with e^ = (— l)^'i+^2+^3+a;4_ inserting these 
restricted fields into @ and dropping the bar on x gives 

b b 

The action (H) with \E' and \E' defined as in (^), reproduces the kinetic term of the action for 
reduced ('real' or 'Majorana-like') staggered fermions, 

Sk = -^Yl VpxXxXx+p ■ (6) 

xp 

The restriction of x and x to odd and even sites corresponds to the projections \E' ^ ^(\E' — 
75^75) and \1' — > ^(\1' + 75^1/75). More explicitly this implies that the matrix fields have the 
structure _ 

The relation of the 2x2 matrix fields ipL,R and ipL r to the fields ip in the target action 
(|I]) becomes clear when writing the Yukawa interaction in (P in the form y{'(pLia'4'R,aj^ij + 
'4'R,ia'4^L,aj4>ij)i where a = 1,2 and i,j = 1,2 are the Weyl spinor and flavor indices, respec- 
tively. 

We can also rewrite the Yukawa term in terms of the matrix fields ^ and ^ if we introduce 
the 4x4 matrix 

/ 



$ = 

I 4)^ 



= -Y.vpip- 



The total fermionic action 

Sf = -T.II1 ^Tr(^.7M*x+A - ^x+a7m*x) + l/Tr(^.v^.$J)] (9) 

X fl 

reduces in the classical continuum limit to eq. (|I]). Using (^ we finally obtain the action in 
terms of the independent x fields 

Sf = -^J2 XxXx+fiiVf^x + ye^CtJix^i^x) = -\Y^ XxM^yXy, (10) 

xfi x,y 

where '^ux — JE ^b ff^ x-b is the average of the scalar field over a lattice hypercube. The 
fermion matrix M in eq. ( ]TI1| ) is antisymmetric and real. 

For completeness we also show the action of the Higgs fields, 

SH = nY: ^Tr(0t <^.^^ + 0t^.0,) - ^ iTr [^t 0, + A(0t 0^ _ i)2] . (n) 

Xfl x 

The total action S = Sh + Sp depends on three coupling constants, the Yukawa couphng y, 
the scalar field hopping parameter k and the quartic self-coupling A. 

We emphasize that this action is invariant under the staggered fermion symmetry group: 
One can check the invariance under shifts, Xx -^ (pxXx+p, ftix ^ (^^^+^(1 - 25^p), 90° 
rotations, axis reversal and global U(l) transformations of the form Xx -^ exp(iae^)xz (cf. 
refs. [0, I). 

The action is not invariant, however, under the full 0(4) fiavor group, and one expects 
to need counterterms to recover this invariance in the scaling region. In the scaling region 
operators with dimension larger than four become irrelevant. There are two operators with 
dimension four which respect the discrete symmetries but break 0(4): 

o« = E^^' o(^) = ;^(^,.+,-yp,.)i (12) 

Xfl X/l 

In general one has to add these operators as counterterms to the action, S ^ S + ciO^^-* + 
€20^"^^ and tune Ci^2 as a function of k and y in order to recover the 0(4) invariance in the 
scaling region. 

Here and also in the numerical work we will restrict ourselves to the case of radially 
frozen Higgs fields corresponding to gSU(2) or equivalently A = +oo. 

3. Phase diagram and fermion mass 

To investigate the phase diagram we carry out a simple mean field calculation in the saddle 
point formulation [|l^] with the replacement v^^x — ^ //^ + ^xf^- The constant fields /^, /**, 
account for a possible ferromagnetic and antiferromagnetic ordering of the scalar field. In 
this approximation we find ^nx ~^ ffi which implies that the fermion fields in eq. (p!OD do 
not couple to the staggered mode /^*. 

At ?/ = the model reduces to the 0(4) model which has three phases: a broken (or fer- 
romagnetic (FM)) phase for k > k^. {k^ = 0.30411(1)), a symmetric (or paramagnetic (PM)) 
phase for —Kc < k < Hc and an antiferromagnetic (AM) phase for k < —Kc- Appropriate or- 
der parameters for a distinction of the various phases are the magnetization t>^ = (^ J2x V^mx) 



Figure 1: The phase diagram at X = cx3. The squares represent the transition points 
determined on an 8'^ lattice, the dashed lines are the results of the mean field calculation. 

and the staggered magnetization f^* = {yJ2x^x^^ix) with magnitudes v = {Y^ii^fiV^'^ and 
v^^ = iJlfii^^ y^'^- In the mean field approximation v^ = /^, f^* = /^*. 

Since the fermions do not couple to /^* the phase transition between the PM and AM 
phases comes out independent of y, fi;^*^~^*^(?/) = —1/4. For the position of the FM-PM 
phase transition we find 

^FM-PM^^ = - - —y' T 5^i^^ ^-- O.OUNDy' , (13) 

where the J2p is a normalized sum over lattice momenta p^ G (— f , f]- We have inserted in 
(|13D the number of SU(2) doublets as the variable N^. In the numerical simulation with the 
HMCA, Nf) has to be chosen as a multiple of two in order to guarantee a positive Boltzmann 
weight, [±(DetM)i/2]Afz,e5« > 0. The mean field results for K^^-^^'\y) and K^^-^^\y) 
{Nd = 2) are represented by the two dashed curves in fig. 1. The lines intersect at the point 
{K,y) ^ (—0.25,4.67). Thus there is a ferrimagnetic (FI) phase where both Vfj, and vf^ are 
nonzero. 

With the mean field ansatz for (p^ ^ the model describes free fermions and it is straight- 
forward to compute the fermion propagator. After Fourier transforming the x SgIcI in the 
action (p!OD and following the steps outlined in refs. [0, [^ we find (using v^ = f^) 

Sab{p) = (5]e'(P+"^)"M-ie-'(P+"^)^)/r (14) 

x,y 

E/.sm p^ + E^,y^vf,cos^P,^ 

with Pfj_ G (— f , f ] and vr^ denoting the 16 momentum four- vectors with components equal to 
or n. The 16 dimensional gamma and fiavor matrices F^ and S^ are defined in refs. 0, ^. 
For p — > we can read off the fermion mass from eq. ([15|) : nip = yv. This reproduces the 
usual tree level relation between rup and v. 



4. Symmetry breaking terms 

Next we want to estimate the magnitude of the 0(4) symmetry breaking terms which are 
induced by the fermions, using renormahzed perturbation theory. To that end we expand 
(iV£)/2)TrlnM in powers of the scalar field. From the two-point contribution we find a 
contact term {5R/2)J2fj.pfi^fi^R, corresponding to the counterterm 0*^^^. For the coefficient 
we obtain 






E^cos>^ 

p Ef^Sm^P^^ + yREf^vlliCOS^Pf, 



R- onyR]^ Z^ Y- o^v.2^ _l.,2y- .,2 



fsNnvl , (16) 



where the subscript R indicates a renormahzed quantity. We remark that keeping the term 
VrJ^ fi 'vfj^^R cos"^ Pf^ in the fermion loop, gives 0{a^) corrections to fs which are negligible deep 
enough in the scaling region but which we want to include if mp = VrVr is of order one. 

The second term er Yl,^ <^t r appears in the four-point contribution to the effective action. 
We find, 

This leads us to consider the following tree-level effective action, 

Seff = -f d'xll E d,^u,Rd,^^A^ + Sr6,.) + ^ E <r + ¥iT. vIr? + ^k E <r] ■ 

(18) 
As a consequence of the £R-term the shape of the effective potential differs from the usual 
Mexican hat form. The potential now has 16 discrete minima at ipR = (±^, ±|, ±|, ±^)vr 
with 

VR = i-ml/iXR + eR)y/'. (19) 

In the infinite volume limit the scalar field will be frozen to one of these minima. Then we 
can replace VRY^/j.vf^jiCos^ Pfj, in eqs. (|16D and (p!7|) by ("^|'/4) X^^cos^p^ and compute the 
prefactors fs and fe as a function of mp. We find fs ~ —0.0150, —0.0110, —0.0070 and 
fe ^ 0.0054, 0.0050, 0.0043 for nip = 0, 0.3, 0.5. 

In order to estimate the effect of the symmetry breaking terms on masses of the longitudi- 
nal (a) and transversal (tt) modes we compute now the scalar field propagator for the effective 
action in eq. (|^). For this purpose we decompose ip into its longitudinal and transversal 
components, ipR^^ = {vr + (jR)e'f^ + tircJ^, where e^ is chosen to point in the direction of one 
of the 16 minima and e^ are three orthogonal vectors. We shall make the convenient choice 



e 



i(l, -1, -1, 1), e^ = i(-l, 1, -1, 1), e^ = i(-l, -1, 1, 1), e^ = i(l, 1, 1, 1). In this basis 



the inverse scalar propagator has the form 



^f)p'S-^ + 2{XR+eR)vlS"'S^'+2eRvlS-^^^+SR[Y.pleX-P^ 



G-\pr^ = il + ^)p'6-^ + 2{XR+eR)vlS"'S^'+2eRvlS-^^^+SR[Y^pleX-^S-'] . (20) 



Since Sr is small compared to 1 we can treat the off-diagonal part (i. e. the last term in 

as a perturbation. Then we find for the propagators G^ijp) = C^^i^p) and 6*^ = 1 J2j G^^{p) 

the forms 

G.Ap) = ^^^yrV + 0{5l)) , (21) 



where the 6j^ correction is bounded by |5|. for all momenta. The particle masses are given 

by 

2 _ 2{Xn + eR)vl , _ 2eRvl 

^ " 1 + 5nlA ' ""--! + 5^4 • ^''^ 

This shows that the three transversal modes in the scalar spectrum acquire a mass as a 
consequence of the symmetry breaking term in the potential. 

We can get an impression of the size of the symmetry breaking corrections for yji < 3, 
Nn = 2 and mp = 0. Then \6r/A\ < 0.068, |4 < 0.055 and Er < 0.88 where the latter one 
has to be compared with typical values oi Xr = 5 — 10. The values which we find for these 
quantities in the physically relevant region are actually smaller than these bounds. It should 
be kept in mind that we have neglected in the relations (|2^) the usual loop corrections not 
involving the symmetry breaking terms. 

Note that to leading order in Sr the propagators G„{jp) and G.„(j>) in (^) are covariant 
functions of p^ and e^, such that we can rotate our frame of reference and choose for example 



e^ = 4. In the following we shall use the approximation (|2l|) for our analysis of the numerical 
propagator results. 

5. Results of the numerical simulations 

In this section we present our results for the phase diagram and give some preliminary 
results for the renormalized couplings hr and Xr. The results were obtained at A = cxo and 
without the two counterterms. The numerical simulations were performed on L'^T lattices 
with periodic boundary conditions in all directions except for the fermion fields which had 
antiperiodic boundary conditions in the time direction. For the calculation of the fermion 
and scalar propagators we carried out the simulations on L^24 lattices with L = 6, 8, 10, 12 
and 16 and at k = 0. 

The phase diagram is shown in fig. 1. The squares represent the transition points which 
were obtained by scanning systematically the (k, y) coupling parameter space in vertical and 
horizontal directions. For each point in the raster we have typically accumulated a statistics 
of 500 scalar field configurations on an 8^ lattice. We used the rotation technique for the 
calculation of the order parameters v and Vst ||TT| , 0. The numerical results are qualitatively 
well described by the mean field results. The phase diagram in fig. 1 is very similar to the 
phase diagrams of other models with a hypercubic Yukawa coupling [O . 



For the determination of the fermion mass mp we have measured the fermion propa- 
gator Sab{p) given in (|l^ for momenta p = (0,^4) and computed from it the projections 
SoipA) = iTr[r4S'(p4)]/16, Sf^ipi) = Tr[S^S5r5S'(p4)]/16. The components S^ (/i ^ 0) are 
non-invariant under cubic rotations of the scalar field and would vanish when averaging over 
many configurations. To avoid finding zero for all of them we have rotated the scalar fields 
such that v^ = f ^^,4. Then, the numerical values ioi Si, i = 1, . . . ,3 are found to be zero 
within errors and we have fitted the non-vanishing components 5*0 and S4 to a free fermion 
propagator form (cf. (p!5|)) 

e r \ Zp^osinpi ZF,AmF,iCOsp4^ 

'J0IP4J -^ jz 2 \ • 2 \ 2~' '-'4IP4J ^ jz 2 \ • 2 \ 2~' ^^"^z 

(1 - mj,Q) sm p4 + mj,Q (1 - mj,J sm p^ + mj,,^ 

where we allow for two different masses and wave-function renormalization constants in 5*0 
and 5*4. In fig. 2 we have plotted Sq^^p^) sinp^ and 5*47^ (^4) cosp4 for various values of y as 



Figure 2: 5*0(^4)3111^4 and 5*4(^4)005^4 as a function of siia^ p^ for several values of y 
(k, = 0, V = 12'^24j. The error bars are in all cases much smaller than the symbols. The 
straight lines were obtained by fitting 5*0(^4) and 5*4(^4) to the forms ^^BJ. 



a function of sin^p4. The data points fall nicely on straight lines which were obtained by 
fitting Sq and 5*4 to the forms (|^). Within the statistical errors, we find that mpfl ~ '"^f,4- 
We therefore have taken the average, which will be denoted by mp in the following. The 
results for mp are summarized in table 1 for several y values and several lattice sizes. We 
find Zp^Q K. 0.8 and Zp,^ approximately 5% smaller than Zp^^. When rotating v^ in a different 
direction we find that mpfl, Zp^ and mp^4^, Zp^^ may differ by about 7% for y < 4.2. Such 
systematic effects are presumably due to scaling violations. 

In order to monitor finite size effects and extrapolate to infinite volume, we have computed 
mp and v on lattices of different spatial extent L. In a previous work with naive fermions 
both observables were found to obey the relations vl = Voo + ai/L"^ and mp^i = "'^^,00 + 02/-^^ 
with ai, 02 > ]12|. This finite size dependence off and mp is due to the massless Goldstone 



bosons. In our model the three transversal modes acquire a mass m^ > as a consequence of 
the 0(4) symmetry breaking and we expect deviations from the 1/L^ behavior if L becomes 
larger than 0{l/m.^). In fig. 3 we have displayed the 1/L^ dependence of v and mp. The 
plot shows that the 1/L^ dependence is well fulfilled for y = 3.8 though small deviations are 
visible for mp on the largest lattice with L = 16. The deviations become more pronounced 
at y = 4.0, which is consistent with the observed increase of mjr (see fig. 4). The effect of 
the symmetry breaking on v and mp, as measured by the deviations from the straight line 
behavior, appears to be smaller than 15%. The effect on the ratio mp/v is even smaller. 

The masses m„ and rriT^ of the longitudinal and transversal modes were obtained by 
fitting the scalar momentum space propagators Go-(p) = {Y.x,y (^xfyye^^^^"^^) /V and Gt,{p) = 



Figure 3: The results for v and mp as a function of 1/L^ for y = 3.8 and y = 4.0. The 
straight lines are drawn to make deviations from a linear 1/L^ dependence visible. 

I J2j{J2x,y vr^vr^e*^^'^~^^)/V for sufficiently small momenta to free boson propagators, cf. (PT|), 

G-^M ^ Z-'Ml,^ + P') , P^O. (24) 

The quantity p^ = 2 I]^(l — cosp^) is a lattice equivalent of the momentum squared in the 
continuum. Here we have chosen the magnetization to point in the four-direction as for 
the computation of the fermion mass (see also the remarks at the end of sect. 4). When 
plotting the inverse propagator G^]^{p) as a function of p^ we expect to find a linear behavior 
if p"^ -C rajp, as is found in the 0(4) model. We measured the a and vr propagators on a 
12^24 lattice where G~\{p) for the three smallest momenta has an approximately linear 
behavior. We note that because of curvature in G~]^(p), the results of a linear ffi tend to 
overestimate m^^T^ and Z^^^^. In a forthcoming publication [^ we will extend the analysis 
of the scalar propagators by taking into account the one-fermion-loop contribution to the 
self-energy which allows for a good description of the curvature in G~]^{p). 

In fig. 4 we have plotted the fit results for m^ (squares) as a function of y. The mass 
has a minimum close to y = 3.8 which is due to two different effects. The increase towards 
small 7/ is a consequence of the finite lattice size as in the pure 0(4) theory. The mass has to 
approach the scalar mass nis in the PM phase continuously when crossing the FM-PM phase 
boundary region. On a finite lattice ms does not vanish at the FM-PM phase transition 
and consequently m^r has to grow when approaching the FM-PM phase transition. The 
increase at large y is due to the 0(4) symmetry breaking. The numerical data for rriTr may 
be compared with the one- loop results given in eqs. ( p^Tf ) and (|^). The values obtained for 
m,T^ are represented in fig. 4 by the diamonds. For y ^ 4.0 the agreement with the numerical 
results is reasonable. 



Figure 4: Numerical results for m^^ (squares) as a function of y (k = 0, V = V1^2A). The 
diamonds represent the values for m.^ which were obtained by inserting the numerical results 
for mp, yR and Vr into eqs. ([7^ and 



We will present now some preliminary results for the renormalized couplings yR and Xr, 
using the tree level motivated definitions yR = mp/vR and Xr = ^{1 + 5ij/4)(mo-/f/j)^ — Er. 



The definition of Xr is based on eq. (22) and we shall use the one- loop values given in eqs. (16) 
and (|17D for 6r and Er. For the normalization of the scalar fields we use the wave-function 
renormalization constant Z^^ (not Z„ because the a particle is unstable). The renormalized 
field expectation value is then defined by vr = v/y/Z^. We find that Z^ is smaller than in 
the pure 0(4) model by roughly a factor four. In table 1 we give the results for yji, \/2Xr 
and the ratio rUa/vR. The table shows that the numerical values for \/2Xr and the ratio 
fn^/vR, which differ by the Er and Sr corrections, agree within 5%. This indicates that the 
effect of the symmetry breaking is relatively small. 

The bare Yukawa couplings are relatively large and it is therefore interesting to compare 
the values for y^ listed in table 1 with the tree level unitarity bound which for this model is 
given by y'^''' = 2Jtx /Nd ^ 2.51 [jrT|. At the edge of the scaling region with rrifj ^ 0.75 the 
numerical results for yji are quite close to this value. This indicates that the renormalized 
couplings are relatively weak, but they appear to be stronger than in the model with naive 



fermions where the numerical results for yR are significantly smaller than y^ |]T2 



The values for the quantity \/2Xr with L = 12 may be compared with the infinite volume 



results obtained previously in the 0(4) model [0. There the ratio J2Xr^oo = 'ma,oo/vR,oo 
ranges from 2.5 to 3.1 for mo-.oo = 0.4 — 1.0. The finite volume results in table 1 with 
rricr^L < 1 lie clearly above this range which indicates that also the upper bound of the Higgs 
mass grows when the Yukawa coupling is turned on. 

6. Conclusion 

We have introduced a fermion-Higgs model with reduced staggered fermions in which the 
staggered fiavors are coupled to the Higgs field. In a numerical simulation with the hybrid 
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Table 1: Some results for nip, v, Z^,, ttIt^, m„, yji, m^/vR and {2XrY/'^ . 
fixed to 0. The lattice size is L^T with T = 24. 



The value of k is 



Monte Carlo algorithm the model contains two isospin doublets in the scaling region. To 
recover the full 0(4) symmetry we have to add two counterterms to the scalar part of the 
action. The addition of counterterms can perhaps be avoided using an F4 lattice. 

In this paper we have studied the model without these counterterms. We find that a 
one-loop computation of the effect of the symmetry breaking can account reasonably well 
for the measured values of m^ and predicts a 5% correction in the relation between m„lvR 
and V2A/J. Also the finite volume dependence of v and mp indicates that the effect of the 
symmetry breaking is relatively small. 

As a preliminary result we find that at k = and at relatively large bare Yukawa 
coupling y the renormalized Yukawa coupling cannot be significantly larger than the tree 
level unitarity bound in the region with m„ < 0.7, rap < 0.5. The quantity a/2Xr comes out 
to be larger than the infinite volume results in the pure 0(4) model. 

The model has proven to be very efficient in a numerical simulation and it is possible 
to perform calculations on relatively large lattices. The time for an update of a scalar field 



configuration is smaller by roughly a factor 10 than in a model with naive fermions [TT], |T2 



We have investigated here a simple example in a class of models in which the staggered 
fiavors are coupled to scalar or gauge degrees of freedom [^, || . We will investigate in future 
more complicated models which might provide a non-perturbative formulation of a chiral 
gauge theory on the lattice. 
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